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A principal motivation to develop graphene for future devices has been its promise for quantum spintronics. 
Hyperfine and spin-orbit interactions are expected to be negligible in single-layer graphene. Spin transport ex- 
periments, on the other hand, show that graphene's spin relaxation is orders of magnitude faster than predicted. 
We present a quantum interference measurement that disentangles sources of magnetic and non-magnetic deco- 
herence in graphene. Magnetic defects are shown to be the primary cause of spin relaxation, while spin-orbit 
interaction is undetectably small. 



Spin lifetimes in graphene are remarkable in that they are 
unremarkable. The first measurements of electron spin re- 
laxation found 100 picosecond lifetimes [1 |-similar to what 
one might expect for conventional metals or semiconductors- 
and more recent measurements confirm that initial result||2][3). 
These experimental values can be contrasted with much more 
favourable theoretical predictions: micro- or even millisec- 
onds are expected[41, due to carbon's low atomic number 
(weak spin-orbit interaction) and the lack of nuclear spin in 
the predominant isotope (weak hyperfine interaction). These 
material properties make graphene an extremely promising 
material for classical or quantum spintronics[5], but potential 
applications await an understanding of the practical mecha- 
nisms of spin relaxation in the material. 

The orders-of-magnitude disconnect between spin relax- 
ation measurements and theory remains one of the most im- 
portant puzzles in graphene research. Recent theoretical work 
has focused on finding a mechanism that would give rise to 
unexpectedly strong spin-orbit interactions, perhaps associ- 
ated with random electric fields from the substrate or localized 
electric fields near adatoms or vacancies. |4| Another possibil- 
ity is that the conduction electrons in graphene are strongly 
coupled to magnetic moments, but this explanation lacks a 
microscopic picture of where these moments arise. It is pos- 
sible to create defects with a magnetic moment intentionally 
in graphene, for example by ion bombardment, [6] [7) but it 
remains an open question whether paramagnetic defects that 
may be present in intrinsic (natural) graphene explain the de- 
phasing rate that is observed at low temperatures. 

Quantum interference is a powerful tool for studying charge 
and spin interactions of conduction electrons with their ma- 
terial host. Random lattice strains, trigonal warping, and 
atomic-scale disorder in graphene dephase the valley sym- 
metry, and dynamic electron-electron interactions give rise 
to inelastic charge dephasing.[8-13| Recently it has become 
clear that additional interactions lead to an apparent satu- 
ration of dephasing at temperatures below a kelvin. lfl2l fl4l 
The precise mechanism has yet to be determined, but two of 
the proposed mechanisms have a direct bearing on the mys- 
tery of graphene's fast spin relaxation: uni-axial spin-orbit 
interactions,! 15 1 or defects with a magnetic moment lfl4l [T6l . 

This paper presents a quantum interference measurement 
that distinguishes magnetic and non-magnetic dephasing 
mechanisms in graphene for the first time. We demonstrate 
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FIG. 1. (a) Simplified device geometry and measurement setup (see 
supplement for device and measurement details). Purple regions are 
graphene, gold indicates leads, and light-shaded regions have been 
etched away, (b) Conductance as a function of gate voltage; the 
shaded region indicates the studied interval, (c) A typical UCF trace 
in perpendicular field, (d) A typical autocorrelation function (solid 
curve) and its derivative (dotted curve, no vertical scale plotted). The 
vertical spike in the dotted curve corresponds to the noise peak in the 
autocorrelation. The dashed line indicates the inflection point, used 
as a measure of coherence through Eq. |TJ(. 



that magnetic moments are indeed a significant source of or- 
bital dephasing in graphene, and are the dominant mechanism 
for spin relaxation|[T]-[3]; in addition, we find a non-magnetic 
dephasing mechanism of similar strength, whose microscopic 
origin is yet to be determined. 

Parameters associated with quantum interference can be 
extracted from a transport measurement in two ways. Uni- 
versal conductance fluctuations (UCF) are seen when the 
phase coherence length is comparable to device dimensions 
(Fig. [TJa,b,c)), and their statistics may be compared to theo- 
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retical predictions. Alternatively, the average low-field mag- 
netoconductance may be fit to weak localization (WL) theory. 
From an experimental point of view WL measurements are an 
easier way to extract a dephasing rate, and historically much 
more common, but it is difficult or impossible to distinguish 
different dephasing mechanisms solely from WL. This exper- 
iment combines UCF and WL measurements to separate and 
quantify the various sources of dephasing in graphene. We 
begin by discussing UCF. 

A monolayer graphene device was exfoliated onto an 
SiC<2/Si wafer (Fig. [TJa)) and measured in a dilution refrig- 
erator with a two-axis magnet. Conductance G was measured 
in a four-terminal configuration, as a function of gate voltage 
VbGi in-plane and out-of-plane magnetic fields Bu and B±, 
and temperature T (Fig. [TJb,c)). Similar results were found 
in a second cooldown (different cryostat) of the same device 
after an annealing step. Conductance fluctuation data for a 
narrow range of Vbg were analyzed by their autocorrelation 
in perpendicular magnetic field, f(5B) (Fig. |TJd)). 

UCF are dephased by any degrees of freedom in a conduc- 
tion electron's environment that change faster than the mea- 
surement bandwidth (hertz). Such dephasing can be charac- 
terized by a rate t^ cf that is the sum of conduction electron 
scattering rates from uncontrollable dynamic sources. fiTl At 
low temperatures, the dominant contributions to t^q F are 
other conduction electrons and dynamic defects in the device, 
whether magnetic or non-magnetic. From an experimental 
point of view, the inflection point of f(8B) provides a robust 
metric of this rate ! 181 : 
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where D = 0.03 m 2 /s is the diffusion constant that was cal- 
culated from G. 

As seen in Fig. 2, the temperature dependence of t^q F 
is linear over the range T = 0.1 • • ■ 1.4 K with a slope 
that is close to the value predicted for electron-electron 
interactions |9] Q3D, and an extrapolated offset t^q F (T=Q) 
that implies a low-temperature saturation of the dephasing rate 
(the corresponding length ^Dt\jcf(T=0) = 1.7/im is much 
smaller than the flake dimensions [Fig. [if a)]). Whereas cer- 
tain symmetry-breaking static environments (e.g. out-of-plane 
spin-orbit coupling) can induce saturation in WL. lfl4l[T5l they 
do not lead to a saturation in UCF. [17 18] Instead, the finite 
t^q F (T=0) observed in this experiment indicates the pres- 
ence of dynamic degenerate defects, that is, defects that do 
not freeze into a single state as temperature is decreased. The 
remainder of this work probes the nature of these degenerate 
defects: are they magnetic, and how strongly do they interact 
with the conduction electrons (how fast do they change state)? 

Performing the UCF measurement with an in-plane mag- 
netic field shows clearly that some of the defects are magnetic: 
dephasing is reduced at B\\ = 6 T (Fig. 121, indicating that the 
magnetic moments have been polarizedto a static configura- 
tion and no longer contribute to t^q F . Quantitatively, the net 
change in with large Bii is the magnetic scattering rate, 
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FIG. 2. Dependence of the rate t^cf on temperature for Bu = (o) 
and Bii = 6 T (•). Dotted lines have slope 9.0 ns /K and offsets 
r -^ p (T=0) = 10.9 ns" 1 (upper line), T ^ F (T=Q) = 6.2 ns' 1 
(middle line). The lower, dashed line shows how t^cf would appear 
without saturation (t^q F (T=0) = 0). 
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FIG. 3. Dependence of t ucf on rms total magnetic field Btot = 
(S| + 2 ' f° r var i ous T. Curves show the theoretical crossover 

for a scattering rate t~* = 5 ns" 1 from spin-| magnetic defects 
with g — 2 (solid), or g = 1 (dotted). 
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This rate is, in itself, an important 
finding, as it corresponds to the spin-flip rate for conduction 
electrons due to unpolarized magnetic defects.] 17] The data in 
Fig. [2] thus prove that magnetic defects induce sufficient spin 
relaxation to explain previous spin transport measurements in 
monolayer graphene. |j}[3j 

The crossover to full defect polarization was analyzed in a 
second cooldown of this device (Fig. [3). The field required 
to turn off the dephasing grows with temperature, as expected 
for the thermodynamics of free magnetic moments. We ob- 
tain the theoretical curves in Fig.|3jby applying the definition 
of Typp in Eq. ([T| to numerically simulated UCF with spin- 1 
defects (see supplement). At high temperatures the behaviour 
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is consistent with defects of g = 2. A significant departure 
is seen at 200 mK and below, indicating the need for a more 
careful theoretical treatment of the magnetic defects as quan- 
tum obj ects . J20l [2T| 

The fact that t^^ f (T=0) does not go to zero at 
high field indicates an additional saturation mechanism 
that is apparently non-magnetic, with dephasing rate 
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^ucf( t =0^||=6 T) a 6(4) ns" 1 [Figs. 
first(second) cooldown]. The data presented so 
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far do not al- 
low us to say more about this non-magnetic mechanism. Is it 
merely device noise that limits UCF? Is it a more fundamental 
inelastic mechanism, such as the two-channel Kondo dephas- 
ing that was predicted for metals a decade ago?[22| Along 
the same lines, it is difficult to ascertain from the UCF data 
whether the magnetic dephasing results from a Kondo-type 
interaction of a few defect spins strongly coupled to the elec- 
tron gas, or from a large number of slowly fluctuating mag- 
netic moments. 

To address these questions we compare the UCF results to 
an analogous measurement based on WL, which is sensitive 
to time reversal symmetry (TRS). Like UCF, WL may be de- 
phased by a dynamic environment, but only when the fluctua- 
tions occur faster than the dephasing timescale — a cutoff time 
nine orders of magnitude shorter than the analogous timescale 
for UCF. Unlike UCF, WL is also dephased by a static envi- 
ronment if it does not preserve time-reversal symmetry; ex- 
amples of this are magnetic fields, or spin-flip processes from 
unpolarized magnetic moments. 

Graphene's magnetoconductance (Fig. |4j a )) is typically fit 
to a WL theory[8| that includes only non-magnetic dephas- 
ing mechanisms, but the UCF data demonstrate that graphene 
also suffers from significant magnetic dephasing. If the mag- 
netic defects vary slowly, they distinguish the spin-singlet 
and -triplet channels of the WL correction. ifTTl l23l l24l which 
complicates the fitting. The dephasing can be more reli- 
ably characterized by extracting the zero-field magnetocon- 
ductance curvature to obtain a single rate Twi, defined as 
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where ^ 



1.05 is the device aspect ratio and G is the av- 
erage conductance. For slow unpolarized magnetic defects 
(B = 0),j23l one expects 
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where T m£ [ g is defined the same as for UCF and r T ^ s is the 



summed dephasing rate from other scattering mechanisms 
that break time reversal symmetry. For fast magnetic de- 
fects, on the other hand, is simply the sum of rates 

As in the case of UCF, is seen to increase with tem- 
perature, and a zero-temperature offset is clearly observed 
(Fig. |4|b)). The common slope with respect to T reflects 




(b) 

















-0.4 0.4 

6 ± (mT) 



0.4 0.8 1.2 
Temperature T (K) 







(c) 














15 














in 












^% 






10 








7 


















6 


~$ 


° a ■ 






5 








5 























0.2 0.4 







Second cooldown, 0.08 K 













1 


2 2 

S n cn 


2 


3 



FIG. 4. (a) Conductance, averaged over the same Vbg range as 
used in Fig. 2. (b) Comparison of characteristic rates r^, L [square, 
extracted from (a)] and r^^ F (o, from Fig.|2| at B\\ = 0. The dotted 
lines are WL and UCF rates expected for electron-electron parameter 
9.0 ns _1 /K, collision rate t~ = 4.7 ns _1 with slow magnetic de- 
fects, and non-magnetic offsets of 3.5 ns^ 1 (WL), 6.2 ns _1 (UCF). 
(c) Dependence of on -Bjf (same Vbg range as used in Fig. pl. 
The solid line is rf^ s (B|| ) = 4.3 ns" 1 + (4.0 ns _1 /T 2 )^, a fit to 
the B\\ > 1 T data. Inset: The same data at low field, plotted linearly 
in B» . 



the equal effect of electron-electron interactions on WL and 
UCF |fl9l , thereby confirming that Eqs. ([T| and Q are not mis- 
calibrated. The effect of an in-plane field on WL (Fig. |4]c)) 
is more complicated than the analogous measurement for 
UCF (Fig. [3} due to graphene's ripples, which convert the 
uniform in-plane field to a random vector potential. This 
breaks time reversal symmetry, ||25l |26l giving T-rrts^ll) = 
t^ s (0)+ /3B 2 , where t^^O) is the inelastic dephasing rate 
from non-magnetic sources and j3 describes the ripple geome- 
try. In the data tZ™ increases sharply from to 50 mT, which 
can be explained by the suppression of two WL channels by 
Zeeman splitting and the resultant transition from Eq. ([3]l to 
t wl = t trs + ! T ma g -GQ) Above 50 mT, t W l decreases at 
first as the magnetic defects polarize and their dephasing ef- 
fect vanishes. For much higher fields (B\\ > 0.5 T) the defects 
are fully polarized and r^,^ nas collapsed to r^g, giving the 
B 2 dependence seen Fig.|4|c). 

Taking UCF and WL data together, we can draw several 
conclusions about the mechanisms of spin relaxation and low 
temperature dephasing in graphene: 
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1. Scattering from magnetic defects induces a spin flip rate | ( 2011) 



. The smaller field-induced sup- 
Fig. 4^c)) is consistent with the 



weaker contribution of r m „, to t wt for slow magnetic defects 



T mag = ^ ns -1 . This is seen directly as the field-induced 
suppression of (Fig. 
pression of T^ h (~ 2 ns" 

l 

mag LU 'WL 

i.e., those that change slowly on the dephasing timescale but 
fast enough to dephase UCF 

2. Spin-orbit interactions can be excluded as a significant 
contribution to spin relaxation in graphene. Spin-orbit cou- 
pling would generate either antilocalization at B\\ = or 
a much larger decrease in t^l for small B», depending on 
the spin-orbit symmetry |[T5l : moreover the temperature de- 
pendence t^q F (T) would be nonlinear] 18 1. None of these 
effects are observed. 

3. WL and UCF each indicate a non-magnetic compo- 
nent to the saturation in dephasing. For the second cooldown, 
the non-magnetic rate for WL dephasing was Tt^ s (T=0) ~ 
3.5 ± 1 ns -1 after subtracting the contribution from electron- 
electron interactions, identical to the corresponding rate for 
UCF, r U( i F (r=0,S||=6 T) ps 3.8 ± 0.2 ns" 1 . This indi- 
cates that the non-magnetic source breaks time-reversal sym- 
metry, and must therefore change rapidly on the dephasing 
timescale. For the first cooldown, the UCF rate was higher 

( t uCf( T= °! B W =6 T ) ~ 6 ' 2 ± °- 3 ns_1) while the WL rate 
was t^ s (T=0) fa 3.5 ± 1 ns -1 ; the higher UCF rate may 
be attributed to defects with dynamics too slow to break time 
reversal symmetry. 

Both magnetic and non-magnetic dephasing mechanisms 
limit coherence in graphene below IK. The magnetic scat- 
tering rate is too large to be explained by remote magnetic 
moments, requiring instead that the magnetic defects are elec- 
tronically coupled to the graphene. Recent WL data|14| sug- 
gests that the magnetic defects may be midgap states at the 
Dirac point, formed at vacancies or edges. For the non- 
magnetic dephasing, we can rule out bistable charge systems 
in the SiC>2 substrate; their broadly-distributed level splittings 
would produce a rate proportional to r. ll22l The data instead 
suggest a class of nearly degenerate non-magnetic defects in 
the graphene itself, whose microscopic origin is yet to be de- 
termined. 
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SUPPLEMENTARY INFORMATION 

M. B. Lundeberg, R. Yang, J. Renard, J. A. Folk 
November 8, 2012 

This supplement covers auxiliary topics relevant to the pri- 
mary text titled Defect-mediated spin relaxation and dephas- 
ing in graphene. The figures in this supplement should be 
viewed in colour. These topics are covered in the following 
sections: 

1. The sample fabrication details are described, as well as 
the measurement and analysis methods. 

2. The expected influences on the device's temperature are 
discussed, in particular the overheating caused by the 
applied bias current. Two measured effects show the 
ability to reach low temperatures: the electron-electron 
interaction correction to conductivity, and the energy 
correlation width of UCF. 

3. We discuss the accuracy of the UCF autocorrelation in- 
flection point when used as a measure of the decoher- 
ence rate. 

4. We describe the classical magnetic defect polariza- 
tion model of UCF, which was used to interpret the 
crossover of t^ cf in Bu . 



Materials and methods 

Single-layer graphene was deposited onto an SiCVSi wafer 
by exfoliation from a natural graphite crystal, then contacted 
by CrAu metallization and plasma etched into a channel. Fig- 
ure [STT] shows the device as imaged by a scanning probe mi- 
croscope. The device was designed for four-probe measure- 
ments, with a channel large enough to be within the quasi-2D 
regime of phase coherence, yet small enough that the conduc- 
tance fluctuations were not self-averaged. In order to avoid 
non-local coherent effects from the metals, the voltage probes 
were etched into long strips of graphene (longer than the co- 
herence length). Note that some damage to the voltage probes 
occurred before the measurements described in the main text, 
so a second CrAu deposition step was necessary. 

After fabrication the device was annealed for 30 minutes 
at 400 °C in N2:H2. The device was measured in dilution 



refrigerators equipped with two-axis magnets, enabling inde- 
pendent control of magnetic field components in the plane of 
the graphene (Bu) and perpendicular to it (B±). Differen- 
tial conductance G was measured by standard four-terminal 
lock-in techniques at 870 Hz (Primary Fig. 1A). The carrier 
density was controlled by the back gate voltage Vbg through 
275 nm of S1O2, giving a typical conductance curve with the 
minimum near Vbg ~ 50 V (Primary Fig. IB). In the first 
cooldown we examined a narrow range of gate voltage Vbg = 
-2.5 ■■■2.5 V (Vbg = -1- • • 1 V for the Bu -dependence 



Current leads Broken probes 




Repair Voltage probes 



FIG. S:l. Phase image from atomic force microscopy, taken after 
deposition of gold leads, plasma etching, and a preliminary measure- 
ment. Electrostatic shock damage is visible as gaps in the graphene 
leads. After this image was taken, the dashed region was filled with 
gold, repairing this voltage probe and allowing the primary measure- 
ments to take place. The channel width and the spacing between 
voltage probes are indicated. 

data), where the carrier density was n s « —4 x 10 12 cm~ 2 
and the diffusion constant was D w 0.030 m 2 /s (estimated 
from the conductivity and density of states, via the Einstein 
relation). A second cooldown of the same device, after a sec- 
ond annealing procedure in N2:H2, produced similar results, 
with parameters Vbg == — 5 ■ ■ ■ 5 V, n 8 ~ —3 x 10 12 cm~ 2 , 
D w 0.033 m 2 /s. The data in the main text and this sup- 
plement are from the first cooldown, except where otherwise 
noted. 

Conductance fluctuation data were analyzed by their auto- 
correlation in perpendicular magnetic field, f(5B), defined as 
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f(6B) = (SG(B ± ,V BG ) SG{B ± + SB, V BG )) 
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B 2 -8B 
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(S4) 
(S5) 



Here, 6G(B±, Vbg) is the fluctuating part of conductance, 
obtained from the raw data G(B±, Vbg) by subtracting off a 
smooth background function. The conductance was scanned 
over a range B\ ~ 50 mT to Bi ~ 150 mT, for ten to 
twenty different values of Vbg spread over the narrow range 
of gate voltage noted above; the averaging over Vbg was done 
to gather more fluctuations and thereby improve the statistical 
accuracy of ( |S5| >. The smooth background, used to obtain SG, 
was obtained by fitting G(B±, Vbg) to a polynomial of the 
form 

go + g\Bx + QiB\ + g 3 V B G + 94 V£g + 05-Bj_Vbg- (S6) 
for free parameters {gt}. 



Device temperature 

There are many claims of saturating physics at low tem- 
peratures, and it is always crucial to establish that the device 
temperature is in fact as low as the temperature of the cryostat. 
Otherwise it may simply be the temperature, not the physics, 
that is saturating. This is of great concern for experiments 
below about 50 mK where it soon becomes difficult to prop- 
erly heat-sink the device in the presence of virtual heat leaks, 
paramagnetic or nuclear heating by a magnetic field, or Joule 
heating due to measurement. 

Regardless of the ease of reaching a 100 mK device temper- 
ature, we have taken several precautions and performed tests 
to remove all doubt about the accuracy of the temperature val- 
ues. Since we examine only temperatures above ~ 100 mK 
in the primary text, and the saturation is already noticeable 
above ~ 500 mK, the following discussion is expected to be 
of supplementary interest. 



Expected influences on device temperature 

The measurements in the first cooldown were carried out in 
the same dilution refrigerator and magnet system as our pre- 
vious experiments of coherent effects in graphene, described 
in Refs. I25landl27l Since those experiments, we have added 
a new wiring and shielding system to the refrigerator. The 
same system has been added to a second cryostat, used for 
the second cooldown. Heat-sinking of the device is ensured 
by providing a 10 mil wiring path from the chip carrier to a 
metal-epoxy-metal heat sink of ~ 200 mm 2 area per wire; this 
10 mil path is free of solder (a superconductor), guaranteeing 



r 



efficient heat sinking at all magnetic fields. High-frequency 
signals originating from the warm parts of the refrigerator are 
removed on every wire by a series of three surface-mount RC 
low-pass circuits (1 .6 MHz, 1 .6 MHz, 16 MHz), followed by a 
2.5 m long lossy RF transmission line; all filters are thermally 
anchored to the mixing chamber so they do not add thermal 
noise. The entire arrangement (filters, device, and heat sink) 
is enclosed by a coaxial radiation shield also attached to the 
mixing chamber. 

In fact, the greatest source of heating in this experiment was 
intentional: In order to expedite the measurement of conduc- 
tance fluctuations with low noise levels, we applied nonnegli- 
gible bias currents through the graphene. The bias current in- 
duces Joule heating, and due to the finite thermal conductivity 
of the graphene this overheats the middle of the device, which 
is the measured region. We calculate the overheating by as- 
suming the device is thermally a one-dimensional system and 
that all heat must leave through the wires by electronic heat 
conduction. Based on the Wiedemann-Franz law, one obtains 
the rms temperature at the hottest point (the center of the de- 
vice), as| 16 1 



V3e 
2nk' 



T — * /^cryostat + ^bias' 



where R is the total electrical resistance from the current 
source's heat sink to the current drain's heat sink, and /bias 
is the rms current. For the Vbg range which was the focus of 
this work, R w 1750 O. 



Figure S:2 shows the variance of UCF for different bias 
currents and temperatures. So long as kT/H is not much 
less than the decoherence rate, the variance is sensitive to 
temperature ifTTl l28ll and so the accuracy of ( |S7| i can be eas- 
ily detected by observing the variance. The figure shows 
how the temperature T « 110 mK (which has UCF variance 
0.43 e 4 /h 2 ) can be produced either by a warm T cryost at ~ 
100 mK or a high bias current Tbias ~ HO mK, as predicted 
by ( |S7| >. The variance saturates below ~ 100 mK, which is 
consistent with kT/h falling below the saturated decoherence 
rate.JTS | Subsection[]explores another aspect of this data. 

In the first cooldown, equation ( |S7| > had its great- 
est relevance for the lowest temperatures T = 
{110 mK,220 mK,310 mK} of the UCF B x -correlation 
data set discussed in the primary text; these tempera- 
tures were achieved by applying /bias = 20 nA with 
Tcryostat = {13 mK, 190 mK, 290 mK}. At higher temper- 
atures we used higher currents (up to 45 nA). Similar bias 
currents were used for the WL measurements (except for the 
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FIG. S:2. Bias overheating effect: the variance, computed from 
G(Vbg) over Vbg = — 5 • ■ • 5 V, as a function of temperature. Here 
B\\ = 0, B ± ~ 100 mT, and 7 bias = 5 nA (o, •) or 20 nA (A, A). 
Open symbols show the uncorrected T cryos tat, and filled symbols the 
corrected temperature T via jS7|> with 7? = 1750 Q. 



110 mK data, which used 10 nA). Similar bias currents were 
used in the second cooldown. 



Test of non-saturating temperature: e-e interaction correction to 
conductivity 

There are other aspects of the device which demonstrate 
that the temperature can be reduced appropriately. One com- 
mon transport test of low temperatures is to examine the 
temperature dependence of ensemble-averaged conductivity 
at nonzero magnetic field. A non-saturating contribution to 
conductivity, proportional to In 77, is expected from electron- 
electron interactions: 



°0 + —r I 1 + c 
irn 



ln(l + F) 



In 77, 



(S8) 



where oo is an offset, F = —0.1 is predicted for graphene 
on SiC>2, and c = 3 at the low temperatures we investi- 
gate here. ll29l [30l This effect should not be confused with 
the electron-electron decoherence which indirectly causes a 
temperature-dependent conductance, via the WL effect. In 
fact, measurements of ( |S8| l should be carried out in a nonzero 
magnetic field to break time reversal symmetry and thus sup- 
press the coherence-sensitive part of WL; this allows a simple 
comparison with dS8). 



In Figure S:3 we show the conductivity over a small range 
in B± near zero, averaged over Vbg - this is a wider range 
of the WL data displayed in the primary text (from the first 
cooldown). Above 1 mT the WL contribution to conductiv- 
ity is insensitive to the coherence and thus is temperature- 
independent; here, WL merely causes a B-dependence of 
the offset Co in Eq. (jS8). Thus, by averaging conductivity 
over a consistent range of magnetic fields above 1 mT, we 
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FIG. S:3. Electron-electron interaction contribution to conductivity. 
(a) VsG-averaged magnetoconductance data expanded from the pri- 
mary text. Shaded regions indicate the field-averaging regions used 
for panel b. (b) Conductivity averaged over both Vbg and B± (cir- 
cles) plotted against log temperature. The dashed line shows the 
slope predicted by theory \S8). 



can test dS8). This VbGi 7?^-averaged conductivity (Fig. S:3 1 



shows the expected logarithmic temperature dependence and 
the slope is close to the expected value. 



Test of non-saturating temperature: UCF thermometry 

As a second test of temperature, we note the recent 
report[31| that UCF themselves can provide information on 
the electron temperature in back-gated devices such as ours. 
In this case, we measure conductance as a function of back- 
gate voltage (not magnetic field), and subtract a linear back- 
ground to yield the conductance fluctuations 5<?(Vbg)' Then, 
Vbg is converted to Fermi energy Ep using the ideal graphene 
relation E p = hvyjira] Vbg — Vq\, where: v = 10 6 m/s is 
the Fermi speed, a — 8 x 10 10 cm~ 2 /V is the gate capaci- 
tance density, and Vq = 45 V is the presumed location of the 
Dirac point in this device. This gives us the UCF in energy, 
Sg(Ep), from which we calculate the autocorrelation function 
in energy: 

f E (SE) = 1 f 2 Sg(E F )Sg(E F + SE)dEp. (S9) 

When kT/h greatly exceeds the decoherence rate t^ 1 , the 
correlation function's shape is entirely determined by thermal 
smearing, so the half-width of f(5E) is expected[31| to be 
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Figure S:5 shows autocorrelations taken at even lower tem- 
peratures, T pa 26 mK according to Eq. at various mag- 
netic fields. The Ei/ 2 correspond to temperatures T w 50 mK 
for the Bii = 6 T data (based on observed high-field dephas- 
ing 6 ns -1 ), or T « 65 mK for the Bu = data (based on 
observed low-field dephasing 11 ns" 1 ). Similar values are 
obtained from E\p. It is strange that the observed widths of 
f{5E) are too large at this low temperature, however we can 
be confident that the temperature is not saturating in the 0.1- 
1 .4 K region examined in the primary text. 



FIG. S:4. 
100 mT. 



Energy autocorrelations of UCF for Bu = 0, B± ~ 
The black curve corresponds to /bias = 5 nA with 



T, 



cryostat 



= 100 mK. The red curve corresponds to /bias = 20 nA 
with Tcryostat = 13 mK. The small peak on top of the black curve is 
an artifact of measurement noise. Correlations were computed from 
G(Vbg) measured over Vbg = —5 ■ • ■ 5 V. [These are two of the 
data points in Fig.|S:2|] 
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FIG. S:5. Energy correlations of UCF for (Bit , B±) values of (0, 0), 
(0, 50 mT), and (6 T, 50 mT). Correlations were computed from 
G(Vbg) measured over Vbg = —5 • ■ • 5 V, with /bias = 4 nA and 



' cryostat 



13 mK. 



El/2 — 2.72fcT, and the inflection point is expected lfl8l to 
be E\p — 2.14fcT. In fact, at low temperatures this condition 
is not satisfied: we have kT/h ~ r^ 1 due to the saturating 
coherence which is the topic of the primary text, so the de- 
coherence rate also influences the width of f(SE). In this 
situation, E 1 / 2 (or E\p) can still be converted to T though in 
a more complicated manner, by taking into account the effect 
of t^ 1 (which is known from -Bip).|[T8) 

Figure [ST4| shows the energy autocorrelation from two con- 
ductance traces in the first cooldown at T ps 100 mK, where 
t^ 1 = 12 ns -1 . By taking into account the value t^ 1 — 
12 ns -1 that was measured in this situation, we compute the 
temperature from the energy widths. [18] In one case the tem- 
perature is mainly due to Tbi as , and its energy widths corre- 
spond to 120 mK (whereas |S7) predicts 113 mK). In the other 
case the temperature is mainly due to T cryosta t, and its energy 
width correponds to 1 15 mK (whereas ( |S7[ > predicts 104 mK). 
The use of either £1/2 or E\p yielded the same temperature 
to within 4%, for these correlation functions. 



Correspondence between inflection point and decoherence rate 



The rate T U( i F defined from the inflection point B\p, 



r^ F = 2eDB IP /3h, 



(S10) 



only corresponds to the true UCF decoherence rate r7 under 
certain conditions, as explained in Ref. 18 We list several 
conditions below; some conditions are nearly violated and 
hence there are some small systematic errors in Def. ( |S 10[ > 
(error meaning the difference between t^ 1 and t^^ f ). Some 
of these errors are known well enough to be corrected, how- 
ever we have chosen not to correct them in order to retain the 
simplicity of definition ( |S10| l. In any case, none of these error 
sources are able to artificially generate the saturation that is 
the topic of the primary text. 



Limiting case: quasi-2D, thermally-smeared 



Definition S10 assumes the quasi-2D thermally smeared 
regime, where the coherence length y^Dr^ is smaller than 
device dimensions, yet is larger than the 'thermal length' 

^/Dh/kT. 

The transition from quasi-2D behaviour! 1 8 1 
( b tp = ^/[ eDT 4>]) to quasi-lD behaviourim 
(B$ = VGh/leWyflh^} for width W) should occur 



in the vicinity of 



D/W . In the present device 



W — 4.3 /im, so the crossover would be when Tj, 



drops 

-1 



below 1.6 ns _1 . The lowest t^ cf measured was 5 ns 
at high Bh and the lowest temperature; this is still a factor 
of three higher than the crossover value, but it is perhaps 
possible that W begins to influence t^q F at this point. We 
are unaware of any predictions of behaviour in the crossover 
and so there is no simple way to quantify this error. As an 
unlikely worst case the crossover could be a gradual quadratic 

crossover (B IP = Bfp ■ [1 + ^Dt^/W 2 ] 1 / 2 ), meaning that 

.-1 



all of our extracted t U( ; f values would be 2 ns higher than 
the actual 1 . 

The thermal smearing limit holds up to 1181 B\p < kT/ eD. 
For the lowest temperatures (~ 100 mK) at zero By, the val- 
ues become as large as B\p = 1.3kT/eD; in this situation 
the error in t\jcf due to the thermal smearing approximation 
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is 10%. For all other data points we have Bjp < 0.7kT/eD 
giving errors less than 4%. 

Background subtraction biases 

The systematic errors originating from the conductance 
background estimation are very small for the inflection point, 
compared to other aspects of the correlation function such 
as variance or half-width. [ 18 1 Based on the considerations 
in Ref. [18] the resulting bias errors in t^ gf are negligible, 
< 3% for the high-temperature data and < 1% for the low- 
temperature data. 

Perturbation by symmetry-breaking disorder 

When there are static symmetry-breaking forms of disorder, 
some (but not all) modes of UCF may be suppressed. If the 
symmetry-breaking rate is comparable to the decoherence rate 
rT , then the correlation function's shape is distorted Ifl8ll32l 
and Def. ( |S 10[ > loses some accuracy. The symmetry-breaking 
rate causes no change if it is much smaller than the decoher- 
ence rate; conversely if the rate is very large then the sup- 
pressed modes are too weak to contribute. 

Graphene's two valley symmetry breaking rates[10| were 
measured by the usual WL fitting formula|8| at fields above 1 
ml We extract r^ 1 = 70 ± 20 ns" 1 and r" 1 = 4000 ± 
200 ns™ 1 , much like seen in other studies. ||9] [25) These 
rates are very large compared to t^q F and hence the valley- 
related modes have little influence. The influence of the r,- 
suppressed mode on t^^ f can be calculated] 18 1 and is less 
than 10% for the high temperature (1.4K) data, and less than 
less than 4% for the low temperature data. The larger rate r" 1 
gives even smaller effects that matter only above 10 K. 

Graphene's spin-orbit rates are expected to be very low 
<C 1 ns -1 and as yet have not been experimentally deter- 
mined. If stronger than expectations, then spin-orbit interac- 
tions would dephase some UCF modes and cause some devi- 
ations in t uQp- A very large (> 5 ns -1 ) spin-orbit rate would 
however produce a specific temperature-dependent shift lfl"8l 
in Tygp which is not observed in the data. 

In conclusion, t^q F measures the rate t^ 1 with typical ac- 
curacy better than 10% for the B\\ = and B\\ = 6 T data. 
There is however another form of symmetry-breaking disor- 
der that is present when Bn is applied, which does lead to 
significant changes in the crossover of 7"uQ F (-Btot); this is the 
topic of the following section. 

Quasi-2D UCF correlations in an in-plane magnetic field 

A large in-plane magnetic field induces a Zeeman splitting 
Ez = 2/is£?tot °f conductance fluctuations, while simultane- 
ously the magnetic defects gain an average polarization. Il20l 
This causes the UCF field correlation function f(SB) to 



evolve in a nontrivial manner as the i?tot is turned on. What 
determines the inflection point t^ cf m ^is case? 

Model 

To answer this question, we simulate the theoretical corre- 
lation function of UCF including high-field effects, and apply 
Def. ( |S 10[ > to predict the value of t^q F . This simulation re- 
quires us to move beyond the one-parameter correlation of the 
primary text to the more general two-parameter correlation 
function of UCFESl EH l33l 

F(SB, 5E) = SG(Bj_,E) SG(B ± +SB,E + 5E), (Sll) 

where SG(B±, E) is the fluctuating part of conductance at 
perpendicular magnetic field B± and energy E. The ordi- 
nary one-parameter magnetoconductance correlation, Eq. 1 
in the main text, corresponds to a measurement of f(SB) — 
F(SB, 0), i.e. it corresponds to SE = 0, since we only com- 
pute the correlation between conductances taken at the same 
Vbg (same E). 

In the analysis of UCF where the spin degeneracy is 
broken (by e.g. spin-orbit interactions ll32l . frozen mag- 
netic impurities [24 1, polarized impurities lEOl . or Zeeman 
splitting [20 1), the full F correlation function is written in 
terms of up to four distinct modes, and each of these modes 
can be written in terms of the spinless correlation F which 
is given by the generic UCF theory. This spinless correlation 
function depends on a single dephasing rate r^ 1 and temper- 
ature T, and is fully specified as F[tT 1 ,T\(8B, SE). In the 
procedure below we use the quasi-2D spinless F computed by 
the procedure in Ref. 18 since our device is in the quasi-2D 
regime. The results can be extended to the quasi- ID case by 
using the spinless correlation function found in Ref. [3T| (with 
the field-dependence noted in Ref.l34l. 

For the case of classical magnetic defects and no spin-orbit 
interaction 1 20 1, 

F{5B,6E) = 2F[t^ 1 ,T](SB,SE) 

+ F[t u \T]{SB,5E + E z ) 

+ F[t^,T](5B,6E-E z ). (S12) 

The first term is comprised of two identical modes, each indi- 
cating the correlation between two conduction electrons with 
the same spin (up/up and down/down). The latter two terms 
express the correlations between electrons with opposite spins 
(up/down and down/up). 

Because of defect polarization in the magnetic field, the 
rates r^ 1 and r^ 1 differ. The rate t^ 1 , indicating the loss 
of coherence between two electrons which are identical (same 
spin), though passing through the graphene at different times, 
is given by[20| 

r n \B tot ) = r^ 1 + [1 - P(B tot f}T m l g . (S13) 
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This may be interpreted as the decoherence rate proper (ie. the 
dephasing by dynamic environment). Here t^ 1 is the field- 
independent part of the decoherence rate (from e.g. electron- 
electron interactions, non-magnetic dynamic defects) and 
T mag i s tne rate °f collision with polarizable magnetic defects. 
P(Btot) is the average polarization of the magnetic defects, 
due to the total field B tot . 

The correlation between conduction electrons of opposite 
spin (but equal kinetic energy) is dephased by a larger rate EUl 



r f f (B tot ) = r^ 1 + [1 + P(B tot ) 2 ]r- 



(S14) 



Note that is always dephased by magnetic impurities, 
even once they have been fully polarized (P = 1); this ex- 
tra dephasing is however not decoherence since it arises from 
static defects. Rather, the extra dephasing in tZV~ simply in- 
dicates that spin-up electrons consistently scatter differently 
than spin-down (this is analogous to the effect of spin-orbit 
interactions in UCF). 

The Zeeman splitting of conduction electrons is also evi- 
dent in Eq. ( |S12| i: a shift of the opposite-spin UCF modes to 
energies SE = ±Ez, so these modes appear as side-peaks 
in the UCF energy correlation function. In the theoretical 
literature these side-peaks are known as m = ±1 diffuson 
triplets [20 21 1. We have directly observed these side-peaks in 
graphene in Ref. 27, and the present device shows side-peaks 
as well, though weakly (Fig. |S:51>. 

Figure S:6 shows how Eq. ( |S 1 2] > evolves as the field is 
increased. Since it is not easy to see the inflection point 
in F itself, we instead plot its second derivative in field, 
^g2F(5B,5E). Again, we remind that f(SB) only mea- 
sures a cross section along SE = 0, thus the inflec- 
tion point Bjp of f(SB) corresponds to the point where 
dSB 2 0) = 0. We draw attention to a few important 

features of Fig. |S:6| 

• For very low fields we have Ez = and P = (thus 
t u = r n ). Here F{SB,SE) = ±F[r^ ,T](8B,8E) 
has the the spinless correlation function's shape. The 
inflection point here is (trivially) a good measure of de- 
coherence, as Ty^p = t^ 1 = x^ 1 . 

• At intermediate fields where Ez < 8kT, the side-peaks 
are close enough to influence the inflection point. This 
causes to be influenced by both r. 
these fields 



Figure 



S:7 



^ and for 



shows precisely how t ucf can 
deviate from the decoherence rate (t^ 1 ) over a range of 
intermediate fields. 



• When Ez is large enough (Ez > 8kT), the larger t^V 
does not matter as the side-peaks are too far away. The 
inflection point coincides with 3Ht^ 



1 /2eD once again, 



thus we measure the decoherence by r, 



UCF- 



The error t vcf /t^ at intermediate fields is influenced by all 
parameters T, Ez, P, t^ 1 , T~* g , so the only way to find its 
value is the direct numerical simulation as described above. 
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FIG. S:6. Second derivative Stt 

ooB 



F{SB, SE) of the simulated UCF 
correlation function ([ST2J, for a case at 3 10 mK. The vertical dashed 
lines indicate the fields Shr^. /2eD and /2eD. (Top-left): 

Near-zero-field case, where Ez = and P — 0. (Top-right): Inter- 
mediate field case where E z = 2.6fcT and P = 71%. (Bottom-left): 
Intermediate field case where Ez = 6.5fcT and P = 98%. (Bottom- 
right): High field case, where E z = YJkT and P = 100%. 



Choice of polarization function 

We have chosen to use a polarization function of the form 
P(Btot) = tanh(<?/iB£?tot/2fcBT), which corresponds to the 
average magnetization of free spin-| moments, normalized to 
P(oo) = 1. One might ask whether this is appropriate given 
that the model used above assumes classical magnetic defects, 
and is not necessarily applicable to spin-| defects. E0l We 
have compared the classical defect model to the much more 
involved quantum calculation[20| that includes the effects of 
inelastic scattering. At least for the central mode (a quantum 
calculation for the side-peaks is not available at this time| 20]), 
the obtained t^q F is nearly the same for the quantum and 
classical models, even for the spin \ case, provided that the 
P = tanh(g hbB \ ot / '2k bT) polarization function is used. 
(Such a correspondence does not extend to WL.) 



Parameters used in the primary text 



In order to generate the simulated curves of Fig. 3 in the 
primary text, we chose r m ^ . = 5 ns _1 for all temperatures, 
based on the typical difference between low and high field. 



11 




FIG. S:7. Comparison between the rate t^^ f extracted from experi- 
mental 310 mK data (first cooldown) and from corresponding simu- 
lation with g — 1.4, compared to the decoherence rate from the 
same simulation. 



The value of r^ 1 , which contains all non-magnetic decoher- 
ence mechanisms such as electron-electron interactions, is set 
by hand in each case to match the high-field t^ cf . The fol- 
lowing table shows the values that were used in the simula- 
tions to produce the curves drawn in Fig. 3 of the primary 
text: 



T (K) 


Va S ( nS *) 


r" 1 (ns- 1 ) 


0.08 


5 


4.8 


0.2 


5 


6.4 


0.5 


5 


9.2 


1 


5 


14.5 



